Generic examples of PT-symmetric qubit (spin 1/2) Liouvillians 
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We outline two general classes of examples of PT-symmetric quantum Liouvillian dynamics of 
open many qubit systems, namely interacting hard-core bosons (or more general XY Z-type spin 
1/2 systems) with, either (i) pure dephasing noise, or (ii) having solely single particle/spin injec- 
tion/absorption incoherent processes. 



PACS numbers: 03.65.Fd, 03.65.Yz, 05.70.Ln 
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The concept of PT symmetry has been introduced [TJ 
as a mathematical framework for studying non-Hermitian 
operators with symmetric (e.g. real) spectra. In recent 
years, PT symmetry has been studied extensively, both 
theoretically [2|-|6[ and experimentally in the context of 
optics @, H| and electric circuits [§| . Very recently, a for- 
mally analogous concept has been proposed [10( on the 
level of master symmetries of quantum master equations 
plTfl3j (Liouville equations) describing open quantum 
systems. It has been shown that the existence of quan- 
tum PT symmetry generically implies the existence of a 
spontaneous symmetry breaking transition, such that for 
sufficiently weak coupling to the environment all coher- 
ences (off-diagonal matrix elements of the system's den- 
sity matrix in the energy eigenbasis) decay with the same 
(uniform) damping rate. 

The aim of this Brief Report is to provide two general 
classes of practically interesting examples to a recent con- 
struction of quantum Liouvillian PT symmetry [loj . We 
use notation exactly as introduced in Ref. [l(| and refer 
to definitions stated there. 

We start by proving a simple observation which seems 
to be crucial for a construction of general examples of 
PT-symmetric Liouvillian systems: 

Lemma. The Liouvillian flow is PT-symmetric, if the 
following conditions are fulfilled: 

(i) The parity superoperator can be represented as 
Vp = UpW, 

where U,W € 8(71) are two unitary operators sat- 
isfying 



and 



U 2 = W 2 = 1, 



[H, U] = [H, W] = 0. 



(ii) 3 M x M real orthogonal reflection matrix Z 6 
0(M,M), satisfying Z 2 = 1 M , such that 



M 



UL m = — S ' Z m rn / L m ,U, 

m' — 1 
M 

WL m = Z m , m ,Ll,W. 



(hi) {i m ,ij„} = c m l, for some c m G K. 
Proof. Using (iii), T>' can be written as 

®P = E ( 2L ™pLl ~ \{[Ll,L m ],p} 



(1) 



(2) 



(3) 



while in the Hilbert-Schmidt metric 

(X>') f P = E (*LlpL m - \{[Ll,L m ],p} 

m ^ 

Then, using (ii), 
V'Vp = V'{UpW) 

= -Uj2 (zL\ nP L m - ±{[Ll,L m ],p}\ W 

Finally, using (i), we have also 

(iadH)V = P(iadH) = -V(iadH)\ (4) 
and therefore 

t'V = -V(C')\ (5) 
i.e., definition of PT symmetry (Eq. (8) of pjj). □ 

One, quite restricted example of application of PT- 
symmetric quantum Liouvllian dynamics has been pro- 
vided in Ref. [10( . Here we show that applications are in 
fact quite abundant and should not be difficult to realize 
in experimentally accessible situations. 

Example 1. Consider n of spins 1/2, or qubits, de- 



scribed by Pauli matrices a 



x,y,z,± 

> 3 - 



1 ... n. Take an 



arbitrary two-spin Hamiltonian with a longitudinal exter- 
nal field 

j<k j 

and up to M = n local Hermitian 'dephasing' Lindblad 
operators 



j = 1 ... 71, 



(7) 
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The interaction strengths J^' k ' z , magnetic field strengths 
hj , and local dephasing rates jj can be completely arbi- 
trary. Liouvillian dynamics defined with respect to the 
Hamiltonian H and Lindblad operators {Lj} is PT sym- 
metric, according to the Lemma above, if we take 

n 

u = n °?> w = l ( g ) 

3=1 

In case J* k = Jj k , hj = 0, the model represents a general 
interacting system of hard-core bosons with zero (or con- 
stant) chemical potential. Then, the Lindblad channels 
([7]) indeed model the pure dephasing noise. 

Example 2. Consider now a slightly more restriced 
Hamiltonian, again for n of spins 1/2, 

H = E(^W°J + J j>N + J lk° z 3°k), (9) 

3<k 

and up to M = 2n Lindblad operators which represent 
incoherent local spin/particle absorption/injection 

L 2 j-i = ajaf, L 2 j = bjaj, j = 1 . . . n. (10) 

Again, the interaction matrices Jj'^' z and the spin- 
flipping rates aj , bj can be arbitrary. The Lemma can 
now be implemented with 

n n 

U=U-h ^ = IK (11) 

In the XXZ-like case JJj. = Jj fe the model repre- 
sents an open interacting hard-core boson model. For 
nearest-neighbor interaction corresponding to one dimen- 
sion (i.e., a chain), and if incoherent processes are only 
at the ends, a,- = bj = 0, 2 < j < n — 1, the model also 
describes an open version of the so-called t — V model of 
spinless fermions. 



In both examples above, the conditions (i-iii) of the 
Lemma are straightforward to check with the reflection 
matrix being trivial Z = 1m- Besides the necessary con- 
ditions [H, U] = 0, [H, W] = 0, implied by (i), we may or 
may not have [U, V] = 0, hence the number of spins n in 
our examples need not be even. 

Provided the parameters in the Hamiltonian can 
be chosen such that both, the energy spectrum and 
the energy-difference (frequency) spectrum are non- 
degenerate, one should observe the spontaneous PT sym- 
metry breaking transition of Liouvillian decay rates while 
increasing the noise/dissipation strength [l0| . and below 
the transition point all the coherences should decay with 
a uniform rate. It is reasonable to expect that both 
classes of examples should be relevant for experiments 
with ultra-cold atom systems. 

We close this note by pointing out the fact that our ex- 
plicit construction (Lemma above) mend the statement 
on the symmetry of the dissipator-perturbation matrix 
V made in Eq. (14) of Ref. [10j . which in general seems 
to be in-conclusive. The symmetry of V is in fact needed 
in order to establish a non-trivial sponaneous PT sym- 
metry breaking transition (i.e., nonvanishing of 7pt)- 
In order to prove Vj.k = Vk,j, we only need to show 
E m \{^3\ L m\^k)\ 2 = E m \{^ k \L m \il)j)\ 2 . Indeed, due to 
the condition (i) of the Lemma the eigenvectors \ipj) of 
H can be chosen to be simultaneously eigenvectors of W, 
say W\ipj) = ujj\ipj), with uij G {±1}, so we have 

^K^|L m |^ fc )| 2 

m 

= J2\^\ WL ^)\ 2 

m 

(1>ALl>wm(il>3\Ll„W\l>k) 

= E \Q>i\iL>\*k)? = E \^k\L m \^)\ 2 , 

m' m 

since Z is orthogonal, Z T Z = 1m- D 
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